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THE CONVERSION OF GODDARD RANGE AND RANGE 
RATE TRACKING DATA TO METRIC UNITS 


ABSTRACT 

The primary purpose of the tracking data conversion described in this re- 
port is to permit tracking system analysis prior to orbit computation. The 
fundamental tracking data recorded during the track of a space vehicle by a 
ground station seldom directly expresses the physical measurements involved. 
In order to directly interpret measured tracking data in terms of physical 
phenomena, prior to orbit or trajectory computation, a conversion of fundamental 
tracking data to metric form is required. The data conversion capability de- 
scribed in this report was implemented to provide a straightforward means of 
converting GRARR tracking data to metric form (i.e., range, range rate and 
angles) free from any data preprocessing such as mathematical smoothing or 
statistical editing. Such metric data is then used as basic information in post- 
flight tracking system analysis. 


•r 


* 


V 


SUMMARY 


Tracking data at the STADAN Goddard Range and Range Rate (GRARR) 
tracking stations is punched on paper tape in the standard Baudot 5 "level tele- 
type code. It is sent to the Goddard Space Flight Center where it is recorded 
on magnetic tape in a format not directly usable by trajectory computation en- 
gineers or scientists. The "STADAN’ 1 tracking data conversion program con- 
verts the GRARR data to MI<S units of range, range rate and angles . The pro- 
cessed data is recorded on magnetic tape. The metric data is also printed out 
with a summary of the data points processed or rejected. 

Data stretches often contain invalid or superfluous data which is of no direct 
interest. The STADAN tracking data conversion program thus also deletes un- 
wanted data stretches; that is, data falling outside of desired timing and tracking 
observable limits. 

The time tagging in this conversion is ground receive time and the desired 
translation to time at the tracked spacecraft is done in the orbit or trajectory 
computation program. The VHF and S-Band GRARR raw data formats are 
described in detail in reference 1. The overall Mission and Trajectory Analysis 
Division tracking data conversion and analysis scheme, of which the STADAN 
program is an integral part, is described in reference 2. 

A thorough discussion of GRARR system sources of measurement error is 
presented in reference 3. 


THE CONVERSION OF GODDARD RANGE AND RANGE 
RATE TRACKING DATA TO METRIC UNITS 

1.0 INTRODUCTION 

The H STADAN n program is designed to provide a means of converting God- 
dard Range and Range Rate (GRARR) data to metric form without such data 
processing as mathematical smoothing or statistical editing. The program 
converts, edits, and reformats tracking data into a standardized format. 


2.0 PURPOSE 

Tracking data from the ST AD AN network is provided to Goddard users on 
magnetic tape in a format not directly usable by trajectory computation engineers 
or scientists. In addition, it often contains invalid or superfluous data which is of 
no direct interest. The "STADAN" program converts the fundamental or ,? raw ,! 
tracking data to units of km, km per second, and degrees, and deletes unwanted 
data stretches, that is, data exceeding timing and tracking observable limits. 

The M STADAN" program computes azimuth (AZ) and elevation (EL) format 
X- and Y-angles in degrees (both are provided to the user), calculates range in 
kilometers, and converts biased Doppler count to average range rate in kilometers 
per second. The processed data is recorded on magnetic tape in a standardized 
format. The rejected data is printed out with a summary of the data accepted 
or rejected. 


3.0 INPUTS 


Table 1 


STADAN Program Input Tape Format 


Character 

Description 

Format 

1-2 

Station number: 

22 -Tananarive (TAN) 

26- Rosman (ROS) 

2 7 - Santiago (SAN) 

28 - Fairbanks (FBK) 
52 -Carnarvon (CRO) 

12 

1 


1 


Table 1 (continued) 


Character 

Description 



Format 

3-5 

GRARR system satellite numbers: 

S-Band 

OGO-A 454 

OGO-B , . . 649 

OGO-C . . , 581 

OGO-D * 673 

OGO-E 814 

OGO-E 951 

GEOS-A 589 

GEOS-B . . 802 

VHF 

IMP-A 346 

IMP-B 460 

IMP-C 542 

IMP-D 658 

IMP-E 770 

IMP-F 751 

IMP-G 953 

13 

6-8 

Day of year 

13 

9-10 

Hour 

12 

11-12 

Minute 

12 

13-14 

Second 

12 

15-19 

X-Angle 

F5.2 

20-24 

Y- Angle 

F5.2 

25-32 

First range measurement 

F8.8 

33-40 

Second range measurement 

F8.8 

41-48 

Third range measurement 

F8.8 

49-56 

Fourth range measurement 

F8.8 

57-63 

First range rate measurement 

F7.0 

64-70 

Second range rate measurement 

F7.0 

71-77 

Third range rate measurement 

F7.0 

78-84 

Fourth range rate measurement 

F7.0 

85 

Band indicator (Cl); 

0 = S~Band frequency (1799.200 MHz) for 3.2 
j MHz subcarrier 

11 


2 


w 


* 


Table 1 (continued) 


«-jrta>*M*3X«sw r *fc *wce* \mn ■ * w *r ' 


(u* ■ » ,-te»;*s 


EE* »T*fcr »«»w« 


Character 


Description 




1 » S-Band frequency (1800.000 MHz) for 2 .4 MHz 

subcarrier 

2 ~ S-Band frequency (1801.000 MHz) for 1A MHz 

subcarrier 

3 - VHF frequency 

4 - PLL operation 


Format | 


86 


Sampling rate and resolution tone indicator (C2): 



Resolution 

Data Rale 

0 ^ 

20 kHz 

l/sec. 

1 * 

20 kHz 

2/sec, 

2 = 

20 kHz 

4/sec. 

3 - 

20 kHz 

6 /min 

4 « 

100 kHz or 
500 kHz 

l/sec. 

5 = 

100 kHz or 
500 kHz 

2/sec*. 

6 « 

100 kHz or 
500 kHz 

4/sec, 

7 = 

100 kHz or 
500 kHz 

6 /min 


II 


87-88 


Range ambiguity 



3.1 Tape Input 

The program reads one input tape file per program run. The format is 
described in table 1. This format is established from the raw data tracking tape 
by the STADAN pre-processor program. The STADAN pre-processor program 
is described in Appendix A. 


3.2 Card Input 

A data card is used to input the time span of data to be selected for process- 
ing. The card is read by the supervisory driver. 
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NOTE 


The Santiago tracking station still supplies range-rate data that is trans- 
mitter frequency dependent. If Santiago data is to be processed, the 
transmitter frequency, in MHz, must be included on the data card. 


4.0 CALCULATIONS 

The program performs the following checks and calculations for each track 
ing data point (record) read from the input tape. Only the data points that pass 
all checks are included in the program output. 


4.1 Time Check 

Time from an input tape record is converted to forms suitable for process- 
ing and output and is checked for validity, or acceptable limits. Only data points 
occurring within the following time ranges are accepted for further processing: 

1 < Day < 366 
0 £ Hour < 24 
0 ~ Minute < 60 
0 ^ Second < 60 


4,2 Azimuth and Elevation Computation 

The following equations are used to compute azimuth and elevation from X- 
and Y-angles: 


EL arcs in (cos X cos Y) 

AZ = arctan (sin X/cosY) 

- arctan [(sin X cos Y/cos EL)/(sin Y/cos EL)] 


( 1 ) 

( 2 ) 


Azimuth follows the usual convention running from 0° to 360° with 0° corre- 
sponding to true North and 90° to East. Elevation varies from 0° to 90°, 0° cor- 
responding to the horizon. 
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4.3 Range Computation (All Stations Except Santiago) 


The following equation is used to compute range in kilometers; 


* 


Range = (Range 


\ 


Raw 



(i<r 3 ) 


( 3 ) 


where; 


Range the computed range in kilometers, 


Range Raw ® the range measurement from the data record which has been 
corrected on-site for station delay. At VTIF, transponder de 
lay has been subtracted on-site. At S-Band, the ’’STADAN 5 ' 
program subtracts transponder delay. 

RA - the number of 8 Iiz ambiguities which have been added to the 

format time interval reading,, characters 87-88 in Table I, 

C = freespace speed of light of 2.997925(10 8 ) meters/second. 

(Reference 4) 


For data before GRARR system modification, the range ambiguity is obtained 
from a separate integrated trajectory program based on simplified celestial 
mechanics which requires an input starting vector. 

The GRARR system has been modified recently to provide unambiguous 
range in which case the integrated trajectory will no longer be required for 
recent data. 


4.4 Range Rate Computation (Santiago only) 


The range rate or spacecraft velocity relative to the tracking station is 
determined from the Doppler shift of a signal transmitted to the spacecraft and 
returned to the ground station. The range rate reading recorded in the data 
format is the number of cycles of reference frequency occurring during the time 
required to count n N n cycles of two-way Doppler plus bias frequency. This 
number, C Q , may be related to the doppler and the average range rate (r) accord- 
ing to the following expression; 
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Table 2 


Bias and Counter Reference Frequencies and Equations Depicting Relationship 
Between C 0 and Spacecraft Velocity. (Reference 1) 



Rosman, Tananarive, Carnarvan and Alaska 

Santiago 

S-Band 

VHF 

S-Band 

VHF 

PLL* 

Crystal 

C rystal 

Crystal 

Crystal 

Inserted 
Doppler Bias 
Frequency 

K 

f 1 1 
4500 

f, 

• 

f, 

500 kHz 

30 kHz 

3600 

5000 

Doppler Time 
Interval Refer- 
ence Frequency 

f 

r 

Pf t 
225 

f, 

180 

f, 

15 

10 MHz 

10 MHz 

Number of 
Refe rence 
Cycles Counted 

— 

c 

0 

/ JL \ 

/ N \ 

N 1 5 

z 

O 

•H 

— 

10 7 N 

225 \ 

/ 180 \ 

1 

5000 \c+f / 

500.000 + ( _)f t 
\C+f/ 

30.000+ )f t 

\C*'r 1 

1 PV 

\4500 \C+r// 

l 1 

\3600 \C r r // 


*p~ phase-locked transponder multiplication constant 
NOTE: See equation 5 for definition of various constants. 


Table 3 


Goddard Range and Range Rate 
"N" Counter Numbers 


Recording 

Rate 

Carrier 

S-Band 

("N" Counter)* 

S-Band 
Count Output 
(Zero Doppler)** 

VHF 

("N" Counter)* 

VHF 

Count Output 
(Zero Doppler)** 

4 per second 

65,503 

1,310,060 

4,093 

1,364,333 

2 per second 


2,620,140 

8,187 

2,729,000 

1 per second 

229,263 

4,585,260 

14,328 

4,776,000 

6 per minute 

3,133,956 

62,679,120 

182,182 

60,727,333 


*Note: S-Band and VHF "N" counter numbers equal the number of cycles of Doppler plus bias 
counted by the "N" counter of the Range Rate Extractor. 

**Note: R Data Format reading at Zero Doppler. 



where 

r = average range rate, 

C = speed of light (2.097925(1 0 8 ) meters/second) 
f b = inserted bias frequency to permit r sign determination (see 
table 2), 

N = N- counter number (a constant, see table 3), 

C average range rate measurement on the input tape corresponding to the 
number of reference frequency, f r , cycles counted in the time required 
to accumulate "N M Doppler plus bias cycles (a variable), 


7 




















f t = transmitted signal from which Doppler is extracted (input on the data 
card). 


Table 4 


STADAN Program Output Tape Format 


Character 

Description 

Format 

1-2 

Station number 

12 

3-5 ' 

Satellite number 

13 

6-7 

Blank 


8-10 

Day of year 

13 

11-12 

Hour 

12 

13-14 

Minute 

12 

15-23 

Second 

F9.6 

24-30 

X- angle 

F7.3 

31-37 

Y-angle 

F7.3 

38-44 

Azimuth 

F7.3 

45-50 

Elevation 

F6.3 

51-62 

Range 

F12.5 

63-71 

Range rate 

F9.6 

72 

System indicator 

11 


The modified GRARR system (Rosman, Alaska, Tananarive, Carnarvon) 
provides a method of Doppler extraction which removes the dependency of the 
velocity measurement on the transmit frequency by deriving the reference and 
bias frequencies from the same source as the transmitted frequency. The 
expression for C c for spacecraft using the GRARR S-Band crystal transponder 
illustrates this independence of the transmit frequency: 


Therefore: 



/n/180 1 \ 

V C o 3600 / 

/N/180 _ _L\— 
\ C„ 3600/ 


(6) 


( 7 ) 
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The values of bias frequency, f b , and counter reference frequency, f r , are 
listed in table 2 along with the equations showing the relationship between the 
data format range rate format number, C o , and the spacecraft velocity. 


4.5 Time Tagging 


Each of the data types of angles, range, and range rate is associated with a 
particular n time tag” (Reference 2). In the GRARR system there is a raw data 
frame time, T p , followed by four frames of data at a data rate T DR . The ap- 
propriate time tag for subsequent frames is, therefore given by: 


T g = T p + kX DR seconds 


,(8) 


where 

T g = frame time 
T f = basic frame time 
T DR = data rate 
k = 0, 1, 2, 3 

The ground measurement time for completion of the range measurement 





( 9 ) 


where 

T m = range lime tag .(ground measurement time) 

T s = frame lime 

T. w - inserted WWV time delay correction 

T r = time increment corresponding to relative phase difference between 
transmitted and received sidetone frequency. 


9 


The ground measurement time for completion of the range rate measurement 
is given by: 


Trm = T s + AT + T w 


( 10 ) 


where 

T^ m s range rate time tag 
T $ = frame time 

T w = inserted WWV time delay correction 

AT ss measured Doppler count time interval required to accumulate N cycles 
of bias frequency plus cycles of carrier phase change. 

Equations (4) and (6), it should be noted here, are stated in terms of 
spacecraft time, and suitable adjustment of the time tagging must be 
made in the orbit determination program. 

The GRAJRR angles are time tagged at ground receive time. 

That is: 


T a = T f + seconds (11) 

where 

T p - basic frame time 

T w = WWV propagation correction 

T a = angle time tag 

Note that only one angle read-out is available in the raw data for every four 
frames of range and range rate data. 

For tracking data recorded after November 1968, WWV correction has been 
made on-site, i.e., T w = 0. 
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69181 8 0 0 69181 10 0 0 START TIME, STOP TIME 

DLOW - 7121.3 MDAY 5 7121.4 


STATION NUMBER HAS CHANGED 
DATA RECORDED PRIOR TO GRARR SYSTEM .MODIFICATION 


28 

951 

181 

4 

25 

29 

-42.29 

71-29 

• 00834004 

00836904 

28 

951 

181 

4 

25 

33 

-40.46 

7 2-29 

• 008457 28 

. 00848706 

28 

951 

181 

4 

25 

37 

-38-57 

73.29 

-00857758 

00860799 

28 

951 

181 

4 

25 

41 

-36.54 

74-22 

. 00870058 

.00873179 


STATION NUMBER HAS CHANGED 


00839826 

-00842767 

5279766 

5285109 

5290400- 

5295653- 

2 

4 

0 

00851702 

- 00854717 

5300857- 

5306016- 

5311131- 

5316194 

2 

4 

0 

00863866 

-00866953 

5321219 • 

5326195- 

5331126- 

5336015- 

2 

4 

0 

00876314 

■00879468 

5340861 • 

5345657. 

5350414- 

5355323- 

2 

4 

0 


26 10 

181 

938 

15-000000 

-75-495 

58-835 

26 10 

181 

938 

15-012744 

0.0 

0-0 

26 10 

181 

938 

15-200104 

0-0 

o.o 

26110 

181 

938 

16-012704 

0-0 

0-0 

26119 

181 

938 

16.200140 

0-0 

0-0 

26210 

181 

938 

17-012663 

0-0 

0-0 

26210 

181 

938 

17 • 200177 

0.0 

0-0 

26310 

181 

938 

18-012623 

0-0 

0-0 

26310 

181 

938 

18-200215 

0.0 

0-0 


329-651 

7-447 

0 0 

0.0 

2 

0-0 

0*0 

1910-27331 

0-0 

2 

0-0 

0-0 

00 

-6- 067524 

2 

0-0 

0-0 

1904-20851 

0-0 

2 

0-0 

0 - G 

00 

6-0588002 

2 

0-0 

0-0 

1898 ■ 15270 

0 0 

2 

0-0 

0-0 

0-0 

-6-049949 

2 

0-0 

0 ■ 0 

1892. 10589 

0-0 

2 

0-0 

o.o 

0-0 

-6 041041 

2 


TOTAL RECORDS READ = 5 


ACCEPTED 


REJECTED = 4 


Figure 1. Sample Printer Output 


5.0 OUTPUTS 


5.1 Tape Output 

The program records the edited, processed data on an output magnetic tape 
in a standardized data format (see table 4). Tapes in this format may be re- 
corded with any density, track level, recording technique, or blocking format that 
can be defined in the JCL of the IBM 360 Operating System. A file mark is 
placed at the end of the data# 


5.2 Printer Output 


The program provides a printout with each program run. The input data 
card is printed. Also printed are the converted time, all rejected records, and 
a summary of the records read,, accepted, and rejected. Figure 1 shows a 
sample of the printer output. 


6.0 ERROR COMMENTS 

Except for the print of all rejected records and the rejected records count, 
no error comments are made. 


7.0 IMPLEMENTATION STATEMENTS 

The program is designed for use on the IBM 360 Computer, Models 75, 91, 
and 95 at Goddard Space Flight Center. It is written in FORTRAN IV for use 
with the FORTRAN H compiler, release 16. 

The FORTRAN I/O logical unit assignments are given under OPERATING 
INSTRUCTIONS. 


8.0 SUBROUTINES 

The program consists of the STADAN (main program), ABNORM, SETCON, 
SETSAN, STACHK, TMECHK, and TRACE subroutines. Figure 2 describes the 
general program flow. The following paragraphs describe the STADAN program. 
Figures 3 through 9 are the FORTRAN listings. 
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Figure 2. General Program Flow 

8.1 ST A DAN Program 

8.1.1 Purpose 

The STADAN program edits and processes tracking data at VHF and S-Band 
frequencies and records the data on an output magnetic tape in a standardized 
format. 

8.1.2 Restrictions 

The tracking data input tape must be in the format described in table 1. This 
format is established by the STADAN pro -processor program described in 
Appendix A. 

Printer output and tracking data output tape must be available. (Refer to 
OPERATING INSTRUCTIONS for logical unit assignments.) 

The input data card must be available. 
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8.1.3 Method 


The program receives the following variables by way of the input data card: 
The start time of data to be processed - 

LOWYR - Year 
I, DAY - Day 
LHR - Hour 
LM - Minute 
LS - Second 

The stop time of data to be processed ~ 


KYR 

- Year 

KDAY 

- Day 

KHR 

- Hour 

KM 

- Minute 

KS 

- Second 


TXF - Transmitter frequency, in MHz, for any Santiago data to be processed. 

The program reads data tape records formatted as described in table 1. 

Each record is tested for the start and stop time span limits. The day-time seg- 
ment of the record is tested for acceptable limits in a call subroutine ABNOBM 
with the following variables as arguments: 

ID - Day of year 

IH - Hour 

M - Minute 

IS - Second 

IF LAG - Data acceptable/not acceptable flag, returned by subroutine 
ABNOHM. 

Upon return from subroutine ABNOHM, depending on the flag setting, the 
record will either be rejected and another read, or accepted and further 
processed. 

An accepted record’s angle measurements are processed and AZ/EL con- 
versions are made. 

A call to subroutine STACHK is made, with the following arguments, to 
determine if the data was recorded before or after the GRARB system modifi- 
cation date. 
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1ST - Station number 

I YE - Year 

ID - Day 

IGOTO - Before/ after modification date flag. 

Upon return from subroutine STACHK, where the flag IGOTO is set to re- 
ject or accept the record, processing either continues or the rejected record 
and the following message are printed before the next record is read: 

DATA RECORDED PRIOR TO GRARR SYSTEM MODIFICATION 

If the data is acceptable, a call either to subroutine SETSAN, for Santiago 
data only, or to subroutine SETCON, for all other stations, is made to set up 
constants for calculations of range rate. 

As noted in the raw data format, there are four range and four range rate 
measurements contained in a record. Each of these measurements is processed 
and stored in arrays. The time element between two records is split into nine 
equal parts. The data is written on an output magnetic tape as follows: 

(1) one record of time-angle measurements 

(2) eight records alternating four time-range measurements with four 
time-range-rate measurements. 

Processing continues until an end -of -file is detected on the input data tape. 
The data summary is then printed as follows: 

TOTAL RECORDS READ = ACCEPTED = REJECTED ~ 

The total records output should equal nine times the accepted records count. 


8.2 ABNORM SUBROUTINE 
8.2.1 Purpose 


The ABNORM subroutine tests the time element of the input data record for 
normal limitations® 

8.2.2 Restrictions 

TV> tracking da; input tape must be in the format described in table 1. This 
format is estal Wished by the STADAN pre-processor program described in Appen- 
dix A. 
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Figure 3. STADAN Program Automated Description and FORTRAN Listing 
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ST 0000 10 
ST00D030 


!*P_ICIT SEAL*a X A-H*0-Z X 
C S T A D A N 

DIWEN&IOM R (41. ROC 4) .ROFI XlA) .1C121 

DIMENSION TRM(A), TR0M<4). RA N! 4 ) . RANO(A) 

niMFMsn>i nini. data uni. ir.nQECini 

OIMENSIDN ICDL! io> 

DATA c. IYR. 9L • Q • I SYS/ 299792. SDO *66 *0*00 *57 *295779 00 *2 / 
. AT A SW IN . 5H? , SDA Y/60.D 0*3 600* 3 0.86 A 00. 00/ 

DATA DCDN/. 2^77770- 3/ 
data MAS 1 /• MAN 1 • / 

HATA T TQACF / * MA !M 1 / 

DA T A CMUS/- 299792.500/ 

DATA RE FREQ /i *0 7/ 

c S T A D A N DATA CARDS A 

r $ 

2. 1 i vvn on _hhmm YY^nn^»iMM,^ tttt.ttt A_ 

C START AND STOP TIME TO ACCEPT DATA • 

r IF THERE IS ANY SANTIAGO DATA. THE TPAKSMITTER * 

r. FREQUENCY MUST B k INCLUDED ON THE DATA CARO. • 

CA_l TRACE! I TR ACE ) 

TTAI 1 r= n — 

IW» I TE= 0 
KC1 = AA 
KC2=4A 
IPEC = 0 
ISTDLD = 999 

ISTCHKs 999 

C YYDDD MMMM SS 6X YYDDD HM MM SS 

READ! 5. 1 1 1LDW YR ,LDA Y.LHR.LM .LS. KYR .KDAY .KHR .KM. KS .T X F 
WRITE! 6. 1 5 )LD AYR ,LD A Y .LHR ,L M ,LS • K YR . KD A Y ,KH K . KM • KS .T X F 

11 FDRN AT! I2.I3.1X.212.1X.I2.6X.I2.I3.1X&212.1X.I2 .2X. F9.3 ) 

IS FORMAT! IX. I 2. I 3. 1 X.2T 2. 1 X.I 2 *6X .1 2 .1 3 .1 X .21 2 «1X . I2.2X .F8 .3 I 

CAL— TM EC HK ! D LO W . LO W YP .1_D A Y . LHR . LM . LSI 

CA_L TMFC HK ! HOA Y.KYR . KDA Y .KHR ,KM. KS> 

WRITE! 6. 12 ) SLOW. HD AY 

12 FORMAT! 1 H 0 • 7MDL0W = • e 8. 1 • Sx, 7HH0AY = • F8.I I 

I YR = .0 WYR 

ZXF*TXF* 1 .06 

IF! ISTDLD I 9 . 777. 9 

9 CONTINUE 


ST 000080 
STQCOflgfl 


• • S/B CALL *A 


ST 000 1 30 
ST 000 1 SO 


ST0C0400 

• IM»UT/DUTPUT 

• lN»UT/OLTPUT 

• 1N=»UT /OUTPUT 

• N»UT/OUTOUT 
TtS/B CALLTT 

A A S/P CALL ww 

ST000A80 A 1N»UT/0UTPUT 
OT AIn»UT/OUTPUT 


siflcasoa 

ST0O0S10 


Figure 3. Continued 


v 


00 


1000 READ! 14* 50 0* END = 7 77 *E RR =350)1 ST ♦! SAT .10.KH.M.1S.X .V • R. RO. IC • I RA 
500 FORMAT! 12. 2! I 3 ) • 3< I 2 ) • 2< e 5. 2 ) .4fF0. 8> .4 ( FT .0 J .2 ( I 1 > • 12 » 


IT ALt_ Y= I TALL Y4 1 
W "* I TE( 6. 498 I 

1 T t. 7 <vA t. rr> .IH.M.1 S.H.Y.Bf 11. 


t ai?\ -am .pit i . t nn i .pc i 2 i . rpi T l» 


c 2PC>< 4 1. IC( 1 ) • IC< 2 ) . I PA 

c 4 98 FORMAT! 5 X. I >. 2! I X.I 3 ) » 3 ( I X ,1 3) .2 ( 1 X ,F5.^ > .4 (1X.F8 .8 ) »4 (U . F7*7 ), 

r 1 2( 1 X. I 1 )• 1 X. I 2) 

CALL A9ND RM ( ID. IH . M. IS. I F LAG ) 

IF( !'.AG .NE • 0 ) GO TD 1000 

TF = i.D«a6AQ C » IH 4 36C0 A — M 4 60 4 IS 

CALL TMFCHX ( O TE ST. I YR .ID.IH.M.IS) 

I F C OTEST *LE • DLO» .OR. OTEST .GE. HDAY ) GC TC 1002 
I c ( IREC .EO. 0 ) GO TO 49 

I FC 1ST .NE . ISTOLO I GO TO 49 
GO TO 4 9 

3GF I TAI I Y-1TA1 i Y»1 

Iw» I TE= I I TE ♦ 1 
WP I TE( 6. 351 ) 

151 FORMAT! 1 H 3 » 2 X » ' R p A 0 "POQR E NC PUNTF RED * / > 


ST0O0720 
STOOD 730 
7AO 

STOCO 7 70 


ST000830 


1002 


GO TO 100C 
CON T IM l»F 


GO TO 1003 

43 IF! I ST .NE • I STOlO ) WR I TE! 6.609) 

*,C9 FORMAT! IH''. 1CX, 26HSTA TI ON NUMBER HAS CHANGED 
IST0lD= I ST 

IF! X .GT. 0.0 1 x = X 4 .005 

IF! * T- X_ = X - .00 5 

IF! y .GT. 0.0 ) Y * V ♦ .005 

IF! Y «L T . C.O ) Y * Y - .005 

XRr X/D 
YQ = Y/Q 

W=D COS! XR )*OCOS! YR ) 

Fi R - PAQ GlMf Ml 

SA7 = DSIN! X- 14DC0S! YR ) /DC D S! E LR> 

CAZ=OSlN! YR ) /DEO S! F.LR ) 

A 7R = PA T AN 2! SA7.CAZ) 


ST000930 

ST0C1C40 
STOC1 050 
ST00I 06C 
STOC 1C70 


AZ= AZR40 
EL = FLR*0 

TFfAZ. T« Q1 AZ =AZ »36Q« QQO 

CAwL STACHK! I ST. I YR .IP .IGOTO) 


4 INR UT/OUTPUT 

• input/output 


•* S/R CALL •• 
•• S/R CALL •• 


m INRUT /OUTPUT 

* INRUT/ 3 UTFUT 

mm SA CALI ** - 

* INRUT/ULTPUT 

* INP UT/O LTPUT 


• •S/R CALL •• 


Figure 3. Continued 





1FC1GOTO) 1203.4203.300 

1203 I«» ITE* I«P I TE+1 

IF( ItfRITE.GT.l) GO TO 1206 

w ° I TE( 6. 1 204 ) • lN*UT/OUTPUT 

120 4 F O RMAT! ?X » * D AT A R Z2 QRJEJ PR IO R TO G ftAHfl SYSTEM MflCIFlCAT ION*/ I * Im> uTin uTPtiT 

120fi *RITE(6.1205HST,ISAT,!0,IH,M,I S,X,Y,9 ,ROiIC.!R* • INPUT /OUTPUT 

12C5 FORMAT! 5X» 12.2! 1X.I 3) .3( 1X.I2) ,2(1X , s 6.2) s 4(lX.F9.3)t4(lX.Fa.Q). 4 INP UT/O UTPUT 

1 2( IX. 11 ), 1X.I2) 

GO TO 100J 
3C0 CONTINUE 

I. c (Kf.l>FQ> i r.l Il.ihi3.x:?.fn.ifi?u r.n rn 

310 KC1= ICf 1 ) 

KC 2= IC( 2) 

I e ( 1ST .NE . 27) GO TO 4<30 

CALL SETSAN1 IC.XNF.31ASFR .SRATL tlBAhO) 44 S/R CALL. 44 

I ST C ►AC - I ST 
T_ J44 

490 C Aw L SETCO^l IC. XNF. SR4TF >IBAND» 4 4 S/P CALL 44 

ISTCHC=IST 

244 je( I ST .NE • I STCHK ) GO TO 310 
I COL 4 = 1. / SRATE 4 .10-6 

I COL 5 = I«4ND / 5 

QO 480 1 = 1.4 

x 1 = I ' 1 STC02130 

TRM ( I)=XI*SRATE 4R ( I ) 

RDF IXC I 1 = 9 »( 1)4 1 .0- 7 

TR0M( I) = TRN(I iTRDFrxt IJ/2. -R ( I ) 

RA= IRA 

RANt I )=f Rf I f /P- 

IFC IST.EO .27) GO TO 482 
AlAMB=( ( XNF/lflO. ) /ROC I 1 1 —DC ON 
OL AMB= AL AMB4CMUS 
RAN DC I ) = FL4MB/( ALAM342. ) 

GO TO 4 80 

4B2 DUKE=I XNF4P EFPEQ l/RDI 1 ) 

CON = B 14 SF9-0UKE 

RAND( I >=C*CON /( 2 • *1 XF-CON) 

480 CONTINUE ST002210 

C THIS BLOCK STORES AND SORTS THE NINE TIME ®CINTS. 4444*4 

I C 3 D E ( i 1 — = — 1 — STOQPPAQ 

T( 1 ) = 0.0D0 


Figure 3. Continued 


— 

J = 0 

DO 2010 I = 2. 
J a J- ♦ V - 
I CD DEC I ) - 2 

HAT A { f \ = D AM< 

8.2 
i \ 


ST 002300 
ST002310 
ST 002 320 
ST002330 
__£T QC.23AQ 

20 10 

T( I | s TP M C J ) 



ST 002 350 


J = 0 



ST0G2360 


00 2020 I = 3. 

9.2 


ST0C2370 


J = J ♦ 1 



ST OC 2 380 


I CO DE ( 11*3 



ST 002 390 


ntTAi i l = PAKini i \ 


<; r nri 2 in n 

20 20 

T( I ) = TP DM ( J 1 



ST00241 0 


I COL 45 = I COL 44 1 0 ♦ IC0L5 




ICDL <11* 

t Z OLA 5 




I C3L ( 2 ) = 

1 COLAS 




1 COL < 3) = 

1 COLAS 




ini a inn 

*• i*m as 




I COL (51 = 100 

♦ I C OLA 5 




IC3LC6) = 200 

♦ IC0L45 




IC0LC7) = 200 

♦ 1 COLAS 




1C0L ( 8 ) = 300 

♦ IC0L45 




ICDL<9» = 300 

♦ IC0L45 




> a 1 



ST002420 


L - 1 



STOO 2430 

20 30 

TM s T( K ) 



ST 002440 


DO 2040 I = K , 

9 


ST 002450 


IF< Tin .LT. 

TM ) GO TO 2050 


ST 00 2 460 


GO TO 2 04 0 



ST 002470 

20 SO- 

TM = Tf 11 



S.TQQ2 A«Q 


L = I 



ST002490 

20 40 

CO *4 T IMUE 



ST 002 500 


IF( L .EO • K ) 

GO TO 2060 


ST 002 51 0 


* 

f* 

II 

** 

J 

►* 



ST OO 2 52 0 


T<< ) = TM 



ST 002 530 


trn = innFM i 



stop? sin 


ICO Dc( _ 1 = (CODE ( K ) 


S T 00 2 550 


ICODEf < ) * ICO 


_ - ■. 

S T 002 560 


OAT s OATML) 



ST 002 570 


DATA ( _ ) * DA TACK) 


ST 00 2 580 


OAT A( K ) = DAT 



S T 0 0 2 59 0 


ICTl a lOlLH ) 





I COL ( L ) = ICOUK) 


Figure 3. Continued 



to 


IC3LOCI = ICO 
2060 K = K ! 

L = K 

I F( < • GE • 9 ) C,0 TO 2070 

GO. 13 -2-03-0 


20 70 CON T I M UC 

00 20 90 I = 1*9 
TEM = T( I ) ♦ TF 

10= TEM / SDA v 

IH = T ( I ) / SHR 

Till = T! I ) - IHY53HR 
M = T( I ) / SM IN 

Till = TC I I - M * SM I N 

1 = icooet 1 1 

1 SAT = ICO! t 1 1 

GO TO ( 301 C. 3020. 3030 > . L 

3010 W» I Tfc ! 1 6. 3600) IST.I SAT. ID .H .m,T(I ) .X.Y. AZ.EL.PL. BL.ISYS 
3SC0 F3BM AT! I 2. I 3. 2X , I 3, I 2 .1 2 . e 9.6 .3 ! e 7. 31 .^6. 3 • F I 2 • 5 . FR .6 • I 1 ) 

C30 10 WRITE! 6. 3601 ) I ST. I SA T • I D .1 H • M .T! I ) • X • Y • AZ .EL • BL • BL. ISYS 
GO TO 2 090 

23020 WRITE! 6. 3601 1 I ST . f SA T ,1 D • I H . m , T ! I > . =3L . B L . BL • BL • C AT A ! I ) . RL . IS Y S 

C36C 1 FORMAT! 3X. I 2. I3.2X.I 3 .1 2.12 . c 1 0. 6 .3 IF 8 • 3 ) .F7.3.F13.5.F10 .6. II ) 

GO TO ? 090 

30 30 WPITE! 16. 3 60C1 I ST . I SA T .1 D . I H . M . T! I I .t3L • BL . BL • BL • BL . OAT A ! I 1. ISYS 

23030 WRITrrf 6. 3601 ) IST.I SAT *ID.IH.M.T!I ) .PL.BL.PL.PL. BL .OAT A! I). ISYS 

20 00 Cl M TIM HE 

IREC = IR EC ♦ 1 
GO TO 1000 
777 COM T IN < IE 

WR ITE! 6. 3 603) I TALLY.IREC .1 WRI TE 

3603 FORMAT! 1H0. 1 CX, • TnTAL RECORDS READ = • . I 5 .1 OX . • ACCEPT EO = ».I5. 

110X. 'REJECTED = » . 1 52 1 

REWIND 14 
WRITE! 1 6. 701 ) 

70 1 FORMAT! 2H- 9 ) 

END FILE 16 
REWIND 16 

STOP 

END 


ST 002600 
ST 002610 
ST 00 2 620 
STQQ263Q 

ST 002640 
ST 002650 
ST0C2660 
STC02670 
SX0Q26a0 
ST0C2690 
ST 002 700 
ST002710 
ST0C2720 
ST002730 


ST002740 


ST002790 


ST002F20 


ST QQ2BS0 
ST 002 860 
ST0C2880 
STC02890 


STQQ291Q 
ST 002920 


• IN 3 UT/3UTPUT 

• inrut/oltput 


y in 3 ut/n uiplll - 


• IN 3 UT/3UTPUT 


♦ inrut/output 

* inrut/output 


• input/output 

• INRUT/OLTPUT 

• INPUT /OLTP UT 

• INRUT/OUTPUT 

• INPUT/OUTPUT 


Figure 3. Continued 



Printer output and tracking data output tapes must be available. (Refer to 
OPERATING INSTRUCTIONS for logical unit assignments.) 

The variables passed as arguments by the STADAN program must be 
available. 

8.2.3 Method 

The subroutine receives the following variables as arguments from the 
calling routine: 

ID - Day of the year, loaded by STADAN 

1H - Hour, loaded by STADAN 

M - Minute, loaded by STADAN 

IS - Seconds 

I FLAG - Good/bad data flag, returned by ABNORM 

Tests are made of the time elements to see if they fall within the normal 
ranges of the individual elements. 

When one is not within the normal limits, the rejection flag is set, the 
record is printed out, and control is returned to the calling program. 

8.2.4 Usage 
— 

The call to ABNORM is: 

CALL ABNORM (ID, IH, M, IS, IFLAG) 

y :i I I! II I 1 1 II 1 1 

8.3 SETCON SUBROUTINE 

i 

8.3.1 Purpose 

The SETCON subroutine determines the constants to be used in calculations 
of range rate for all stations except Santiago. 

8.3.2 Restrictions 

The tracking data input tape must be in the format described in table 1. 

This format is established by the STADAN pre-processor program described 
in Appendix A. 

» > V " I 
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SU0PR3 GRAM A9N0RM 


PROGRAM TYPE - SUBROUTINE 


INPUT/OUT 3 UT 

TYPE OPERATION FORTRAN LOGICAL UNIT TYPE ACCESS 

WRITE «... ■ ft SEQUENTIAL, 


» 

THIS SUBPROGRAM 

HAS 19 SOURCE CARDS. 


ADDITIONAL ENTRY 

POINTS? NONE 





COMMON REQUIRED? 

NONE 



SUBPROGRAMS REQUIRED? TRACE IDCQMW 


-THE TOTAL PROGRAM LENGTH IS 566 BYTES. 


SUBROUTINE ABNORM! ID. IH , M, IS. I FLAG I 00010 

IMPLICIT REAL * 8 ( A~H, O-Z I 00020 

C IFLA G. .= 0 IMPLIES GOOD DATA 00040 

C IFLAG .NE. 0 IMPLES DATA OUT OF NORMAL RANGE 00050 

DATA rWir- /> aamp • / .. 

CALL TRACE! I TRACE ) 

IF;_ AG = 0 00060 

I F ( 10 .IT, 1 .OR. ID .GT. 366 ) GO TO 77 00070 

IF ( IH .LT. 0 .OR. IH .GT. 24 ) GO TO 77 00080 

IF ( M .LT. 0 .OR. M .GT. 60 ) GO TO 77 00090 

LEI — I.S— «L.T-^ 0 — .. I.S-..GJ.,— ftQ,.. I—G-O — T Q-77 , CO 100 

GO TO BO 00110 

77 IFLAG = I 00120 

WRITE! 6, 600 ) ID. IH, M, IS 00 130 

C WRITE! 6. 3601 ) I ST. I. SA T .1 D .IH . M .T ! I ) .X . Y . AZ .EL. BL. BL. IS YS 

C 360 1 FDRMAT! 5X. I 2 . 2X , 2 1 3 ,2 I 2 , F 9. 6 ,3 { F 7. 3 > , F6. 3 . F I 2 . 5 . F9 .6 , I 1) 

. ..ftPO FTRM.A.T.! — l.HC, 1QX. PAH-TIME. .IS OU.T TV : RANGE .1 0-X ..13.3.12— )— „ GO 140. 

00 RETURN 00150 

END 00 160 


** S/P 


CALL ** 


* INPUT/OUTPUT 


* IN? UT/O U.TRUX 


Figure 4. ABNORM Subroutine Automated Description and FORTRAN Listing 
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Printer output and tracking data output tapes must be available. (Refer to 
OPERATOR) INSTRUCTIONS for logical unit assignments.) 

The variables passed as arguments by the STADAN program must be 
available. 

8.3.3 Method 

The subroutine receives the following variables as arguments from the 
calling routine: 

IC - two-dimensioned array loaded by STADAN 

IC(1) - contains the character which denotes either S- Band or VHF. 

IC(2) - contains the character which denotes the sampling rate and the 

highest resolution tone used. 

XNF - a variable containing the S-band or VHF ”N” counter number 
returned by SETCON 

SR ATE - a variable containing the recording or sampling rate, returned 
by SETCON 

IBAND - a variable indicating that the data is S-Band or VHF, returned 
by SETCON. 

Upon entering subroutine SETCON, the values of IC determine the constants 
to be used by program STADAN in calculations of range rate. 

After the constants are set, control is returned to the calling program. 

8.3.4 Usage 

The call to SETCON is: 

CALL SETCON (IC, XNF, SRATE, IBAND) 


8.4 SETSAN SUBROUTINE 
8.4.1 Purpose 


The SETSAN subroutine determines the constants to be used in the calcu- 
lation of range rate for the Santiago tracking station. 
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SUBPROGRAM SETCON 


PROGRAM TYPE - SUBROUTINE 


NO INPUT/3UTPUT 

IS USED IN THIS 

SUBPROGRAM. 



... THIS SUBPROGRAM 

HAS 2C SOURCE 

CARDS. 

ADDITIONAL ENTRY 

points: none 


COMMON required: 

NONE 

— ■ 

SUBPROGRAMS REQUIRED: TRACE 




--THE TOTAL PROGRAM LENGTH IS 7i4 BYTES* 


SUBR3 UT INE SE TCONI 1C * XNF . SP ATE . I BAND) 

IMPLICIT REAL*8( A-H ,0-Z > 

DIMENSION Id 4) .RECORD ( 8) . SC0UNT(4) . VCOUNT (4 ) 

DATA RECORD /I O 0 . . 53 0 . . 25D 0 . 1 0. DO .1 . DO . . 5 DO , . 25 DO . 1 0 . DO / 

DATA SCn UN TZ» 22.92 63D 6 .., .LSI f) 07i~) 6 ..-0655Q3 Dfi .3 . 1 A3QSf.nft/ 

DATA VCOUNT/. 01 43280 6,. 0081 Q7D 6 ,. 004093 06 ,. 182 182 06/ 

DATA ITRACE/*SETC V 

CALL Ta ACE( 1TRACE ) *■* S/P 

I BAN 0=0 
J 1= ICC 1 ) + l 

-IEXJJL».fiE.*A.) IBA.N.D.s.5 _ . 

J 2- IC( 2 ) + 1 
SRATS=RECORD( J2 ) 

I F ( J 2.GF..5) ISR=J2-4 
IFIJ2.LE.4) X SR = J 2 
LNF= ISR + IBAND-I 

-T.E.M .NF...I E-. 4 J . XMF-jzSCnUNTT J SR 1 , „ 

IFlLNF.GE.5) XNF = VC OUNT( I SR) 

RETURN 

END 


Figure 5. SETCON Subroutine Automated Description and FORTRAN Listing 


call** 
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8.4.2 Restrictions 


The tracking data input tape must be in the format described in table 1. 
This format is established by the STADAN pre-processor program described 
in Appendix A. 

Printer output and tracking data output tapes must lx? available. (Refer to 
OPERATING INSTRUCTIONS for logical unit assignments.) 

The variables passed as arguments by the STADAN program must be 
available. 

8.4.3 Method 

The subroutine receives the following variables as arguments from the 
calling routine: 



IC 

- two-dimensioned array loaded by STADAN 


IC(1) 

- contains the character which denotes either S-Band or VHF 

1 

IC<2) 

- contains the character which denotes the sampling rate and the 
highest resolution tone used 


XNF 

- a variable containing the S-Band or VHF ”N M counter number, 
returned by SETSAN 

BIASFR 

SRATE 

- a variable containing the bias frequency, returned by SETSAN 

- a variable containing the recording or sampling rate, returned 
by SETSAN 


IBAND 

- a variable indicating that the data is S-Band or VHF, returned 
by SETSAN. 


Upon entering subroutine SETSAN, the values of IC determine the constants 
to be used by program STADAN in the calculation of range rate from Santiago. 

After the constants are set, control is ieturned to the calling program. 

8.4.4 Usage 

The call to SETSAN is: 

CALL SETSAN (IC,XNF, BIASFR, SRATE, IBAND) 


» 

e 
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SUBPROGRAM SCTS*N 
PROGRAM TYPE - SUBROUTINE 


,—NCL. IMP UT/O U TP U T .IS-. USED- IN THIS SUBPROGRAM., 


— THIS -SUBPROGRAM HAS 26 SOURCE CAROS. 


ADDITIONAL ENTRY POINTS: NONE 


COMMON required: NONE 


SUBPROGRAMS REQUIRED : TRACE 


-THE TOTAL-PROGRAM LENGTH IS 790 BYTES. " ~ 


— — SUBROUTINE. SETSANMC.XNF .BIASFR .SRATE.1RANDL 
IMPLICIT REAL*8< A-H ,0~Z) 

DIMEN SION 1C.I 2) .R SPED ( 8). • SRAD AR (A) . VRADLQ (4 ) .BFUG EC2 ) 

DATA R SPED/1 .DO. .SO 0, .2500 *1 O.DO.l * DO ..5 DO o .2500 . 10 .DO/ 

DATA 1 711 nn?nA.. _ 

DATA VRADIO/. 014328D6.. 0081 87D6..004093D6. .18218206/ 

— DATA ITRACE/ ? SANT ' / 

CALL TRACE! ITRACE ) 44 S/R CALL ** 

... BIASFRs Q.ODO 
I BAN 0= 0 

J U-ICt i l.T-I . 

IFCJ1.GE.4) IB A ND = 5 
J 2= IC! 21+1 
SR A TE=R SPED! J2) 

IF! J 2.GE.5) I SR-J2-4 . 

IF(J2.LE.4) ISR=J2 

LN Es LSaLAIBAN DtaJ ' . . 

IF! LNF .LE .4 ) GO TO 10 . 

1F!LNF.GE.5) GO TO 12 
10 XNF= SRADAR ( I SR ) 

BIASFR= 5000 00. D 0 
GO TO 15 

1 2 - ,.. XN F.S.V 3 A n i nl-I SR .) , 

0 1 A SFR = 300 00.00 
15 RETURN 
END 


Figure 6. SETSAN Subroutine Automated Description and FORTRAN Listing 
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8.5 STACHK SUBROUTINE 


8.5.1 Purpose 

The STACHK subroutine determines whether the data was recorded before 
or after the date of modification of the tracking stations to the GRARR system. 

8.5.2 Restrictions 


The tracking data input tape must be in the format described in table 1. This 
format is established by the STADAN pre-processor program described in 
Appendix A. 

Printer output and tracking data output tapes must be available. (Refer to 
OPERATING INSTRUCTIONS for logical unit assignments.) 

The variables passed as arguments by the STADAN program must be 
available. 

8.5.3 Method 

The subroutine receives the following variables as arguments from the 
calling routine: 

I~T - a variable containing the station number, loaded by STADAN 

IYR - a variable containing the year, loaded by STADAN 

ID - a variable containing the day, loaded by STADAN 

IGOTO - flag set depending on whether the data was recorded before or 

after the modification date 

Upon entering subroutine STaCHK, the modification date for each station is 
determined. The day of the year in the data record, passed as an argument, is 
tested against the particular stations modification date. The record is rejected 
in subroutine STADAN if the data was recorded prior to that stations modifica- 
tion date. (The modification year is set at 1969. All data prior to this will be 
rejected.) 


NOTE 

The Santiago station will not be modified. Therefore, all Santiago data 
will be accepted up to this point of flow of the program. 

After testing is completed, control is returned to the calling program. 
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SUBPROGRAM STACHK 


PROGRAM TYPE - SUBROUTINE 


NO INPUT/3 UTPUT 

IS USED IN THIS SUBPROGRAM. 



THIS SUBPROGRAM 

HAS 16 SOURCE CAROS. 

ADDITIONAL ENTRY 

points: nonf. 

COMMON REOUIREO; 

NONE 

SUBPROGRAMS REOU IREO : TRACE 


THE TOTAL PROGRAM LENGTH IS 55fl BYTES* 


SUBROUTINE STACHKt I ST . I YR , l D . I G OT 0) 

DATA I TRACE/' STCKV 

CALL TRACE! ITR ACE) **5/R 

JYR=69 

IGQ.XQa.-l 

IDC HK= 0 

IF! 1ST .EG. 22) IDCHKslAS 
IF< IST.EO .26) I DCHK = 12 0 
LFl IST.EO .2S) IDCHK”298 
I F < IST.EO .52) IDCHK=2A3 

— IF( I ST . F Q .??.) I CCLTH s] ...... 

IF( I GO TO ) 3. 3. 

A IF( LYR.LT.JYR ) IGOTQ-O 
IP( 10 .L T. IDCHK ) I GO TO =0 
3 RETURN 
END 


Figure 7. STACHK Subroutine Automated Description and FORTRAN Listing 


CALL** 
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8.5.4 Usage 


The call to STACHK is: 

CALL STACHK (1ST, IYR, ID, IGOTO) 

8.6 TMECHK SUBROUTINE 

8.6.1 Purpose 


The TMECHK subroutine converts the time element of each record plus 
the start and stop times of the data to be processed, (supplied by the data 
card) to more readily usable numbers. 

8.6.2 Rest r ictions 

The tracking data input tape must be in the format described in table 1. 
This format is established by the STADAN pre-processor program described 
in Appendix A. 

Printer output and tracking data output tapes must be available. (Refer to 
OPERATING INSTRUCTIONS for logical unit assignments.) 

The variables passed as arguments by the STADAN program must be 
available. 

8.6.3 Method 


The subroutine receives the following variables as arguments from the 
calling routine: 

D - a variable containing the converted time, returned by TMECHK 

IYR - a variable containing the year, loaded by STADAN 

ID - a variable containing the day of the year, loaded by STADAN 

M - a variable containing the minute, loaded by STADAN 

IS - a variable containing the second, loaded by STADAN. 

Upon entering subroutine TMECHK, the time elements from the data records 
as well as the acceptable data time span from the data card are converted to 
more read ;, y usable numbers. This makes the task of testing the data time 
against the time span of data to be piocessed easier to perform. The test is 
done in STADAN program. 
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SUBPROGRAM T¥CCM« 


PRQGAtM Trot - SUBROUTINE 


NO XMPUT/O UTPUT IS USIO IN THIS SUBPROGRAM. 


THIS SUBPROGRAM HAS 17 SOURCE CAROS. 


ADDITIONAL ENTRY POINTS: NONt 


common required: none 


SUBPROGRAMS REQUIRED: TWACE 


Tne TOT At. PROGRAM LENGTH IS S2 0 BYTES. 


SUBROUTINE TMECHKID .1 VR.10.IH.M.IS) 

IMPLICIT REAL •«< A-H ,0-7 » 

DATA X TRACE TMCK •/ 

CALL TRACE! I TRACE ) •• S/R 

r»P» iHPiftnn.M* ada» s 

DFspR/86400. 

NYRS= IYR-50 
«« 0 
J« 52 

IS IF< IYO- J ) 50 , SC • A 0 

AC i'Mt, A| 

J* J 

GO TO 1 5 

SO DO*N YR S* 3 AS AH TfC- 1 
DsDDTO^ 

RETURN 

END 


Figure 8 


"MECH*' Subroutine Automated Description and FORTRAN Listing 


call •• 
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After converting the time element, control is returned to the calling 
program. 

8.6.4 Usage 

The call to TMECHK is as follows: 

For the start-stop time span - 

CALL TMECHK (DLOW, LOWYR, LDAY, LH, LM, LS) 

CALL TMECHK (HDAY, KYR, KDAY, KH, LM, KS) 

For the data record time - 

CALL TMECHK (DTEST, 1YR, IH, M, S) 

8.7 TRACE SUBROUTINE 

8.7.1 Purpose 

The TRACE subroutine is used primarily for debugging purposes by keeping 
tra^ 1 ' of the sequence of subroutine calls. 

8.7.2 Restrictions 

The variables passed as arguments by the calling subroutine must be 
available. 

8.7.3 Method 

The subroutine receives the following variables as an argument from the 
calling routine: 

I A - a variable containing the literal data constant from a data state- 

ment in the calling routine. 

The literal data constants are stored in an array dimensioned 90. When the 
array is completed, that is after 90 calls to subroutines have been made, it is 
dumped on unit 6. For normal runs, the unit is dummied. (Refer to OPERAT- 
ING INSTRUCTIONS.) 
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SUBPROGRAM TRACE 


PROGRAM TVPe - SUBROUTINE 


INPUT/nuT»UT 

type OPERATION »QRTHAN logical unit 

TYPE access 

«R|TK 6 

SEQUENT IAI, 

This SUBPROGRAM mas 19 SOURCE CAROS. 


AOC IT TONAL aNT»V POINTS! NON* 


COMMON REQUIRED ! NUNL 

• • 

SUBPROGRAMS REQUIRED! I 3C UM • 



THE TOTAL PROGRAM LENGTH IS HI2 BYTES. 


SUBROUTINE TRACE ( IA) 

DIMENSION IB1 901 

DATA I/O/. IBLK / • •/ 

1 A< I .EQ. 0 ) GO TO 111 

100 CQMTlMUg 

IB1 I ) « IA 
I * I ♦ 1 

IF! I .GT. 90 ) GO TO 200 
GO TO 7 77 

111 00 112 J « 90 

ill latU-A.uu 

I » 1 

GO TO 100 
200 I * 0 

RPITEC 6,6001 I0( J).J = 1 .901 
600 FTRMATl SX, 3 OA 4 > 

JZZI f-OMTINUF 

RETURN 

ENO 


AM ENO (IF DATA SE T HAS BEEN 2 KC CIUNTC RED • 


7 StlOPWlGWAM? HA VK 0 r f. K PROCtSSzO. 


1 02H4, IS THE TOTAL LENGTH >*- ALL TH t ABCVe MOOULES. 

PROCESSING « ILL NO* T*w*MNATt£. 

Figure 9. TRACE Subroutine Automated Description and FORTRAN Listing 


* in^ut/uutput 

* INPUT /OUTPUT 
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8.7,4 Usage 


The call to TRACE is: 


CALL TRACE (ITRACE) 

9.0 OPERATING INSTRUCTIONS 

9.1 Logical Unit Assignments 

The following FORTRAN logical units are assigned: 

Logical unit 5 - card input 
Logical unit G - printer output 
Logical unit 14 - tracking data input 
Logical unit 1G - tracking data output 

9.2 Program Deck Setup 


The program deck consists oi the source deck, data card, fixed JCL cards 
that vary with the run parameters. A sample program deck is given below. 
Consult references 4 and 5 for a complete description of JCL cards needed for 
all possible run parameters. 

Sample Deck Setup: 


// JOBCARD 
// EXEC FORTRANH 
// SOURCE. SYSIN DD* 

. . . Source Deck . . . 

/* 

// EXEC LINKGO 

// GO. FT 14 F001 DD UNIT = 2400 - 9, DISP = (OLD, KEEP) LABEL = (1, BLP), 
// DCB = (RECFM = FB, LRECL = 88, BLKSIZE = 8800, DEN = 3), 

// VOL = SER = XXXXXX 

// GO. FT 16 F001 DD UNIT = 2400 - 9, DISP = (NEW, KEEP), LABEL = 

u. blp), 

// DCB = (RECFM = FB, LRECL = 72, BLK SIZE = 7200, DEN = 3), 

// VOL = SER = XXXXXX 
// GO. DATA 5 DD * 

. . . Data Card . . . 

/* 
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APPENDIX A 


STADAN PRE-PROCESSOR PROGRAM 


INTRODUCTION 


The STADAN pre-processor program is used to reformat the tracking data 
into an acceptable format for processing by the STADAN program. 

PROBLEM S TATEMENT 

Tracking data from the STADAN network is provided to Goddard users on 
magnetic tape in a format not usable by the STADAN program. See table A-l for 
the teletype transmitted format. 

The STADAN pre-processor program eliminates special characters, deletes 
records not in proper sequence, and reformats the data into one record. 


INPUTS 
Tape Input 

The program reads one input tape file per program run. 
Card Input 

There are no data cards required. 

OUTPUTS 
Tape Output 


The program records the edited, reformatted data on a standard label out- 
put magnetic tape. See table A-2 for the output format. A file mark is placed 
at the end of the data. 

Printer Output 

The program provides a printout with each program run. The first two 
input records ar~ printed plus an output record count. 
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Met hex! 


Contained in the input data records are two special characters which denote 
sequence. The first line of the message should have a period preceding the time 
of year and the second line should have a slash preceding the identification 
data. When this sequence is detected, the special characters are deleted and the 
data is written on the output magnetic tape in the single record format described 
in table A-2. 


O PE RATING I NSTKUC T IONS 


Log ical Unit Assi gnments 

The following FORTRAN logical units are assigned: 

Logical unit 5 - card input 
Logical unit G - printer output 
Logical unit 8 - tracking data input 
Logical unit 9 - reformatted data output 

Program Deck Setup 

The program deck consists of the source deck, fixed JCL cards, and JCL 
cards that may vary with the run parameters. A sample program deck is given 
below. Consult references 1, 2 and 3 for a complete description of utilities 
programs and JCL cards needed for possible run parameters. 

Sample Deck Setup: 


// JOB CARD 

// EXEC PGM = IEHINITT 

// SYSPRINT DD SYSOUT = A 

// LABEL DD DCB - (DEN * 3), UNIT - (2400 - !), 1, DEFER) 

INITT SER = XXXXXX, NUMBTAPE = 1, OWNER = 'XXXXXXXX', DISP 
REWIND 

/* 

// EXEC FORTRANH 
// SOURC E.SYSIN DD* 

. . . Source Deck . . • 

/* 

// EXEC LINKGO 

//GO. FT08F001 DD UNIT = 2400-7, DISP = (OLD, KEEP), LABEL = 
(1, BLP), 
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Table A-l 

STADAN Pre- Processor Program Input Tape Format (GRARR Data Format) 


CO 

oo 


TIME 


LI 

L2 


BIT 


W1 


C/R l/F F/S XXXXX 
C/R L/F F/S YYYYY 


W2 

Q/D RRRR.tRRR 
Q/D RRRRRRRR 


2 3 4 5 C 7 8 . 1 23456789 

- - — 

RANGE 

DATA 


W3 


• RRRilRRR 

Srrrrrrr 


W4 


TIME OF YEAR DATA 


D 10 D . HR u> HR , MIN .o MlN i SEC ..i SKC 1 


/ SAT SAT SAT STA STA Cl C2 


i 12345678 . 1 2 


3 4 5 6 7 8 


RANGE 

RATE 

DATA 


IDENTIFICATION DATA 


Ra 


Ra 


W5 * 

Q/D.RRP.RKRRR 
Q'D raiRRRRRR 


10 i 1 23456789.12345678. 


-we- 


vuurumn 

•RRHRRRR 


RANGE 

DATA 


RANGE 

RATE 

DATA 


C/R CARRIAGE RETURN 
L/F UNE FEED 
F/S FIGURE SHIFT 

X ANTENNA POSITION (SIGN AND FOUR DECIMAL DIGITS) 
Q/D QUALITY DA;' A (SPACE INDICATES ALL LOOPS LOCKED, 
RANGE TONES ON, ANTENNA AUTO TRACK) 

R RANGE RATE (0000000-9°99999) 

SEC SECONDS (00-59) 

MIN MINUTES (00-59) 


R RANGE (00000000-99999999) 

HR ‘HOURS (00-23) 

D DAYS (000-365) 

Y ANTENNA POSITION (SIGN AND FOUR DECIMAL DIGITS) 
SAT SATE LUTE IDENTIFICATION (000-999) 

STA STATION IDENTIFICATION (00-99) 

Cl SPACE CRAFT FREQUENCY 
C2 SAMPLE RATE (0-4) AND RESOLUTION 
Ra RANGE AMBIGUITY • (00-63) 


NOTE: DATA IS PUNCHED IN STANDARD BAUDOT 5- LEVEL CODE. 

THERE ARE 52 CHARACTERS IN EACH UNE, MAKING A TOTAL OF 104 
CHARACTERS TO BE TRANSMITTED SERIALLY. 

•SPACE = ARC NOT IN USE; ? = ARC NOT ACQUIRED 

g 0 - SAMPLE RATE AT 6/MIN; SPACE = SAMPLE RATE AT 4. 2, OR l/SEC 


// DCB * (RECFM = U, LRECL = 50, BLKSIZE « 50, DEN = 1, TRTCH = FT), 
//VOL SER XXXX 

// GO.FT09F001 DD UNIT = 2400-9, DISP = (NEW, KEEP), LABEL = (I, SL), 
// DCB ■ (RECFM 1 I), LRECL ■ 88, BLKSIZE 8800, DEN 3), 

// VOL = SER = XXXXXX, DSN = SDATA 1 

/* 


Table A-2 


STADAN Prc- Processor Program Output Tape Format 


Column 

Description 

Designator 

1-2 

Station number 

STA 

3-5 

Satellite number 

SAT 

6-14 

Time (day, hour, minute, second) 

T 

15-19 

X-angle 

X 

20-24 

Y-angle 

Y 

25-32 

First range measurement 

R1 

33-40 

Second range measurement 

RU 

41-48 

Third range measurement 

R2 

49-56 

Fourth range measurement 

R22 

57-63 

First range rate measurement 

RR1 

64-70 

Second range rate measurement 

RR11 

71-77 

Third range rate measurement 

RR2 

78-84 

Fourth range rate measurement 

RR22 

85 

Band indicator 

Cl 

86 

Sample rate and resolution tone indicator 

C2 

87-88 

Range ambiguity 

RA 









SUBPROGRAM MAIN 


MHIIONAM l*Mt - MAIN 


I NPUT/OUt PUT 

MPi ORFPATION MIMHAN CUBICAL LAN I I TVPF Al. CASS 


■ Ml ft 

M» AO 

•Mire 

«NOf ILF. 


6 stoorsfiAL 

M SkOUINllAL 

V SI OOF NT I AC 

T 


THIS SUHPROGMAM has AX SOUHCt CAPOS. 


IMPLICIT I NTtGEMAXI A-/ I 

O I MISSION LlN e |< All »L INI X I AT I . STA( XI .SAT ( J) . T (Ml »•«(*)( .CIXMlA) • 

I VV| SI. M| ( A) .Mill Ml . -»XI 11 .MXXt 8) .MM | < ?» , MM| | | 7) ,«■»< T| . MM XX I M 
OATA 'IOOTTiSLASH.HIANK/'h •••/ ».* 

C <JU I VAl I NCI (LlNfcllb). Qll.lt INI I I 2 i I . *M I I . I L I NF 1 I J 3 I , 0 I 1 I 

EOUl VALENCE I L I Next A I ,0X1 , | L INC342 J I .»»2 I • «L INf 21 J )> .02 2 I 

F QUl VAl I NO (l INI II t .» . * Ml I I I • IL INF II 7) ,M | | | ) I . IL I Nt I I II. I .MM I | III 

HWI VAl INC! I I I NT I I 24 A • T I | I I , I L I Nf I I JA I . » I I I I ) > . I L I NF I I A 1 I * MM | l | l | ) 

UJUIVALI Nil (L INtXt I I *W| I I I . IL INF XI n.M2lll>.lLlNF2llbl.MM2ll)> 

F QUl V Al F NCI IL I Nf XI XA i . SA T I | ) ) , I L INEXI 27 ) , ST At I ) I 
I OU I VAl I NCI (L INt 21 2 T> .C I 2 MI I I I 

r UU I V A L I NC I I L In* XI <A > .M2X I I > > • I L INE XI AJ I . MM22 I I ) I 
|MfC»0 

6 Ml ADI 1, | O.INO.NO IL I Nl I 
| b FUMMA T I AN A | | 

IF I S«»l .FQ. S*. ASM I GO TO b 
HI SPI ULANk I GO TO fo 

GO TO I A 

| A 'll Ai)( H , | 0 i ' N,)i NO ll INIX 

IM SPX .F i. SI ASm| GO TO XA 
IMSP2.F O. BLANK ) GO TO b 
OO 16 I » I • AM 
lb L Ini I I I I «L INF XI I I 
GO TO I A 

X A IIIQI.N' •IILANK) GO TO b 
|F|0| | .Nt .HLANKI GO TO b 
IF I OX. NF . HLANK I GO TO b 
ir ( oxx .Nt . ML ANA I GO TO b 
IF I IMFC.GT .0 I GO TO l<T 
• P I IK b. X7 I L I NF I 
■ M | TT I b, X 7 | l I NF X 
X 7 FOMMAf | SK , «QA| / I 

IT AR I Tl I T , ISISTA.SAT . T » XX. V V * M | .HI | . M X • R 2 2 t MM | « H Ml 1 . MMX . MP 22 .C I XM 
I S FORMAT | HHA I I 
IREC« |HtCF I 
GO TO b 

TO ■»! TEIb.XbllMEC 

x b format i iox.'aftfm ioiting and consolidating, therf aehe •.is.*heco 

IMDS AMITTFN.*) 

END FILE T 

STOP 

ENO 


• I NPU T/OUTPOT 

• | NPOTXOJTPOT 


• | NPIJT/OUTPUT 


• I NPUT XOOTPOT 

• I NPUT/OOTPUT 

• I nput/output 

• I NPOT XOO TPUT 

• I NPOT/JUTPOT 


• I NPU T/OUTPOT 

• I NPUT/OUTPUT 

• I NPUT/OUTPUT 


F igure A- 1 . 


STADAN Pre-Processor Program Automatic Description and FORTRAN Listing 
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